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An unsteady Reynolds-averaged Navier-Stokes solver for unstructured grids is used to compute the 
rotor airloads on the UH-60A helicopter at high-speed and high thrust conditions. The flow solver is 
coupled to a rotorcraft comprehensive code in order to account for trim and aeroelastic deflections. 
Simulations are performed both with and without the fuselage, and the effects of grid resolution, tem- 
poral resolution and turbulence model are examined. Computed airloads are compared to flight data. 
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Nomenclature 

value at the blade tip 
fuselage angle of attack [°] 
advance ratio 
azimuthal position, [°] 
density [slugs/ft 3 ] 
solidity 

collective pitch, [°] 

negative of lateral cyclic pitch, [°] 

negative of longitudinal cyclic pitch, [° 

speed of sound [ft/s] 

local blade chord [ft] 

thrust coefficient 

rotor thrust [lb] 

Mach number 


c/4 sectional axial (chord) force coefficient ( , ) 

jpa^c 


M z C m sectional pitching moment coefficient ( x M \ 2 ) 


2pa 2 c 2 ' 
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M 2 C n 

sectional normal force coefficient ( T Fn 0 ) 

v 2pa z c J 

M x 

rotor rolling moment [ft-lb] 

My 

rotor pitching moment, [ft-lb] 

R 

rotor radius, [ft] 

r 

radial position, [ft] 

BDF 

Backward Differentiation Formulae 

BL 

baseline mesh 

CAMRAD II Comprehensive Analytical Model of Rotor- 
craft Aerodynamics and Dynamics II 

CFD 

Computational Fluid Dynamics 

CO 

coarse mesh 

CSD 

Computational Structural Dynamics 

HRLES 

1 Hybrid Reynolds -Averaged Navier-Stokes/Farge- 
Eddy Simulation turbulence model 

RB 

refined-blade mesh 

SA 

Spalart-Allmaras turbulence model 

SST 

Shear Stress Transport turbulence model 



Introduction 

The prediction of rotorcraft airloads presents a significant 
challenge for Computational Fluid Dynamics (CFD). The 
combination of complex flow physics, nonlinear structural 
dynamics, and the geometrical complexity of rotorcraft all 
must be accounted for. Significant progress has been made 
in the last decade toward meeting this challenge. Potsdam et 
al. (Ref. 1) combined the structured-grid CFD code OVER- 
FLOW (Ref. 2) with the rotorcraft Computational Structural 
Dynamics (CSD) code CAMRAD II (Ref. 3) in a loosely- 
coupled fashion and applied the coupled system to the UH- 
60A Blackhawk helicopter. A code such as CAMRAD II is 
also referred to as a ‘comprehensive’ code, because in ad- 
dition to structural dynamics models, aerodynamic, propul- 
sion, and trim models are provided. In a coupled CFD/CSD 
system, the CFD solver replaces the aerodynamic analysis 
of the comprehensive code, while the CSD code continues 
to provide the structural dynamics and trim. The generally 
good comparison with flight-test data obtained in Ref. 1 
renewed interest in the loose-coupling approach originally 
proposed by Tung et. al (Ref. 4). Indeed, the loose-coupling 
method has now been applied by many other researchers, 
not only to the UH-60A, but to other configurations as well. 

Most of the CFD/CSD work to date has been performed 
using structured-grid CFD solvers. Structured-grid solvers 
are generally very efficient in terms of computational re- 
sources, but for complex geometries, the human effort re- 
quired to generate the requisite grid systems can be signif- 
icant. Unstructured-grid CFD solvers require greater com- 
putational resources, but the human effort to generate grids 
about complex geometries is typically far less. As a result, 
the total time to solution can be less for unstructured-grid 
solvers. In Reference 5, the unstructured-grid CFD solver 
FUN3D (Ref. 6) was loosely coupled to the CAMRAD II 
code and applied to the HART-II configuration for three 
test conditions. Abras (Ref. 7) coupled the FUN3D solver 
to the DYMORE CSD code and applied the resulting sys- 
tem to UH-60A high-speed and low-speed flight conditions. 
Sankaran et al. (Ref. 8) used a novel unstructured/cartesian 
hybrid solver coupled with the RCAS comprehensive code 
to compute rotor airloads for several UH-60A flight condi- 
tions. 

In this paper, we apply the FUN3D/CAMRAD II anal- 
ysis system to the UH-60A for two flight conditions, 
one high-speed and the other high-thrust. The UH-60A 
database (Refs. 9, 10) has an extensive range of aerody- 
namic and structural-loads data at numerous flight condi- 
tions across the flight envelope and thus provides a very 
useful data set for validation of computational results. The 
focus here is on the validation of computed rotor airloads, 
and the assessment of the sensitivity of the airloads to sev- 
eral CFD-related inputs: grid size, turbulence model, and 
time step. The sensitivities of the off-body aerodynamics 
(e.g. vortex propagation) to these variations is not consid- 
ered. Although the emphasis is on aerodynamics, because 


the aerodynamics and structures do interact with one an- 
other some comparisons of structural loads from flight and 
computation are also presented. However, no attempt is 
made to assess the sensitivity of the results to any structural- 
modeling input. 

CFD/CSD Methodology 

The unstructured-grid flow solver used for this study is 
FUN3D. The code solves the Navier- Stokes equations, in- 
cluding the full viscous terms (i.e. the thin-layer assump- 
tion is not made), with several turbulence modeling options. 
FUN3D is a node-centered code, so that the number of un- 
knowns is directly related to the number of nodes in the 
grid. The code supports ‘mixed-element’ meshes in which 
the nodes may be connected into prismatic, hexahedral, or 
pyramidal elements instead of, or in addition to, tetrahedra. 
The solver has a variety of mesh-motion options, including 
rigid, deforming, and overset meshes, and a robust implicit 
time-advancement scheme (Ref. 11). For overset meshes, 
the DiRTlib (Ref. 12) and SUGGAR++ (Ref. 13) codes are 
used to facilitate communication between disparate zones 
in the mesh. These features make the solver attractive for 
rotorcraft applications, where the flow is fundamentally un- 
steady and the flexible rotor blades undergo large motions 
relative to the fuselage. O’Brien (Ref. 14) first applied the 
FUN3D solver to rotorcraft simulations, with the restriction 
of rigid blades with prescribed motion. Subsequent modifi- 
cations resulted in a more general rotorcraft capability, with 
the ability to account for aeroelastic effects and trim via 
coupling with the CAMRAD II rotorcraft comprehensive 
code (Ref. 5). The coupling is implemented via the loose 
coupling strategy outlined in Ref. 1 , and thus is appropriate 
to steady, level flight. 

The rotorcraft CSD code used for this study is CAM- 
RAD II. The aerodynamics modules within CAMRAD II 
are based on lifting-line models utilizing airfoil tables, cou- 
pled with wake models. Although such aerodynamic mod- 
els can provide reasonable results for many flight condi- 
tions, in some cases the predictions of the airloads can be 
inaccurate because of limitations of the relatively low-order 
aerodynamic modeling. The goal of the loose coupling ap- 
proach is to replace the low-order lifting-line aerodynam- 
ics of the CSD code with the higher-fidelity Navier-Stokes 
aerodynamics of the CFD code. Within CAMRAD II, each 
blade is modeled as a set of nonlinear beam elements. In 
addition to the structural dynamics modeling, CAMRAD II 
offers a sophisticated trim capability. For the UH-60A sim- 
ulations in this paper, a three-degree-of-freedom trim is uti- 
lized, with the (solidity-weighted) thrust coefficient, pitch- 
ing moment and rolling moment specified as trim targets 
within CAMRAD II. 

In the loose coupling approach, airloads data from the 
CFD solver and blade motion data from the CSD solver are 
exchanged at relatively infrequent periodic intervals, for ex- 



ample once per revolution or more generally once per inte- 
ger multiple of the blade passage. In a typical coupled simu- 
lation, the initial execution of FUN3D is carried out for two 
complete rotor revolutions, using blade deformations from 
a trimmed CAMRAD II solution with unmodified lifting- 
line aerodynamics. In subsequent coupling cycles, FUN3D 
is run for two blade passages between coupling cycles (i.e. 
one-half revolution for a four-bladed rotor such as the UH- 
60 A). Although coupling could be performed as frequently 
as each blade passage, running the solver for the extra blade 
passage helps damp out transients between coupling cycles. 
This coupling strategy has not been examined for optimal- 
ity; however, it has proven to be a robust strategy (stable 
and convergent). The coupled system is tied together by 
separate interface codes that perform data translation be- 
tween FUN3D and CAMRAD II. A shell script is used to 
orchestrate the execution of the various codes, and to pro- 
vide restart, post-processing, and archiving functions. 

CFD Grids 

The baseline composite-grid systems for the isolated UH- 
60 rotor were developed by over- setting component grids 
for each individual blade with an empty (background) 
square domain. The component unstructured grids for 
the overset computations were generated with the VGRID 
v4.0 advancing-layer and advancing-front grid generation 
software package (Ref. 15). The grids generated with 
VGRID were fully tetrahedral. However, VGRID uses an 
advancing-layer technique to generate the boundary-layer 
portion of the grid so that prisms can be reconstructed in the 
boundary layer for use with a mixed-element discretization. 
In the boundary layer, three tetrahedral cells are combined 
into one prism. The mixed-element grids have essentially 
the same number of nodes and nodal spacing as the fully- 
tetrahedral grids although the number of cells and the shape 
of the control volumes differ in the boundary layers. 

The blade surface geometry used for grid generation is 
based on an updated definition provided by the Sikorsky 
Aircraft Corporation via the UH-60 Airloads Workshop in 
May 2009. The updated geometry includes more span wise 
resolution at the tip cap and a correction to the spanwise 
location of the trim tab. The near-body blade grid extends 
away from the blade to a cylindrical outer boundary of ra- 
dius 1.5 c t ip. In the wall normal direction, the grid spac- 
ing is set such that an average normalized coordinate y + is 
less than one for the first grid cell at the wall for the ma- 
jority of the blade. The maximum spacing at the blade grid 
outer boundary is approximately 0. 10c r//; . The characteris- 
tic spacings of the baseline rotor grid are summarized in Ta- 
ble 1. Note that these are the spacing values associated with 
the VGrid “sources” which are used to control the unstruc- 
tured grid spacing field, rather than the actual spacings be- 
tween node points, though they are loosely correlated. The 
maximum spanwise stretching is 30:1 near both the blade 
leading edge (LE) and trailing edge (TE). The surface grid 


has approximately 5 points across the blunt trailing edge. 
Each component blade grid has approximately 2.65 million 
nodes. Figures l(a)-(c) show a comparison of the surface 
mesh near the tip of the rotor blade for the coarse, baseline 
and refined blade grids. Figures l(d)-(f) show the corre- 
sponding slices taken through the blade volume grids at the 
r/R = 0.865 radial station which is at the outboard edge of 
the trim tab. 



Direction 

Characteristic Spacing 

Blade 

chordwise 

0.0029c LE 

Surface 


0.0014c TE 


spanwise 

0.0009R Root 


at mid-chord 

0.0056R Mid-Span 
0.0009R Tip 


normal 

2.62 x 10 ~ 6 c tip 
(geometric stretching 1.20) 

Blade 

isotropic 

0.10 c tip 


Outer 

Boundary 


Table 1. Characteristic spacings of baseline blade grid. 

The off-body grid (background grid) for the isolated ro- 
tor is defined in a square box whose sides extend 5 R out 
from the rotor hub. The finest spacing in the off-body 
grid is approximately O.lOc^. This minimum spacing is 
maintained, within the constraints of the unstructured mesh- 
ing software, in a cylindrical volume which extends 1.1R 
in the blade plane and 0.20R above and 0.20R below the 
blade plane. This cylindrical volume of refinement was 
achieved by utilizing VGrid volume sources (Ref. 16). The 
background component grid has approximately 7.0 million 
nodes. The baseline composite grid has a total of approx- 
imately 17.6 million nodes. A slice through the baseline 
isolated UH-60 composite grid in Fig. 3 shows the spacing 
characteristics of the blade and off-body grids through the 
computational domain. The grids are shown in a vertical 
plane passing through the rotor hub location. Note that the 
fringe and hole points are not plotted. 

Although a mathematically consistent grid refinement 
study in the sense of Ref. 17 was not performed, alterna- 
tive composite-grid systems were developed to look at the 
effects of grid refinement: a globally coarsened grid and a 
refined blade grid. A comparison of global characteristics 
for these composite grids is provided in Table 2. The coarse 
composite grid was developed from the baseline grid by in- 
creasing the minimum off-body mesh spacing to 0.20c^ 
and increasing the baseline grid surface spacings by a fac- 
tor of 1.3. This coarse composite grid has approximately 
6.9 million nodes. A refined blade grid was developed by 
globally refining the blade surface-grid spacings by a factor 
0.73. The maximum spacing at the blade grid outer bound- 
ary was set to 0. 10 c t i p so that this grid could be overset with 
the baseline off-body grid. 





Baseline 

Coarse 

Refined 

Blade 

Blade Grid 

Nodes (millions) 

2.7 

1.3 

7.6 


Viscous Surface Nodes (thousands) 

62 

33 

126 


Surface Spacing Factor 

1 

1.3 

0.73 


Wall Normal Spacing Factor 

1 

1.3 

1.00 


Characteristic Outer Boundary Spacing 

0.10 c tip 

0.20c^ 

0.10 c t ip 

Background Grid 

Nodes (millions) 

7.0 

1.8 

7.0 


Finest Spacing 

0.10 ctip 

0.20 c tip 

0.10 cup 


Outer Boundary Extent 

5R 

5R 

5R 

Composite Grid 

Total Nodes (millions) 

17.6 

6.9 

37.4 


Table 2. Comparison of composite grid systems. 


An off-body grid containing the UH-60 fuselage body 
was developed to study the influence of the body aerody- 
namics on the rotor airloads. The fuselage surface geom- 
etry was provided by Sikorsky Aircraft Corporation. The 
rotor hub, shaft, and blade linkages are not modeled, nor is 
the tail rotor modeled. Figure 2 shows the surface mesh for 
the fuselage and the baseline rotor grid. The surface mesh 
for the fuselage body has 59,000 nodes. In the wall normal 
direction, the grid spacing (3.80 x 10 c t i p ) is set such that 
an average normalized coordinate y + is less than one for the 
first grid cell at the wall for the majority of the body. The 
outer boundary of the fuselage grid is a square box whose 
sides extend 5 R out from the rotor hub. This grid also con- 
tains a cylindrical volume of refinement in the vicinity of 
the rotor plane of motion similar to the one defined for the 
isolated rotor. The maximum spacing within this cylinder 
is not allowed to exceed O.lOc^. The fuselage component 
grid has approximately 8.4 million nodes. When the fuse- 
lage grid is combined with the baseline rotor blades, the 
resulting composite grid has a total of approximately 19.2 
million nodes. A slice through the baseline rotor/fuselage 
composite grid in Fig. 3 shows the spacing characteristics 
of the blade and off-body grids through the computational 
domain. The grids are shown in a vertical plane passing 
through the fuselage centerline (the fringe and hole points 
are not plotted). Note that when the fuselage is included, 
the rotor is inclined at the shaft angle, while for the iso- 
lated rotor, the rotor is not inclined, and the shaft angle is 
accounted for in the specified angle of attack. 

For all composite grid systems, SUGGAR++ is used in 
a pre-processing step to create the composite grid and the 
domain connectivity file (DCI) for the nominal rotor po- 
sition. For dynamic overset cases, the DCI file is created 
at each time step (azimuthal location) by calling a set of 
SUGGAR++ library routines directly from FUN3D. SUG- 
GAR++ offers a selection of hole-cutting algorithms. For 
this work, the direct-cut method was used to define the 
blanked points in the composite grid, and an overlap min- 
imization process was used to identify fringe points. A 
donor quality value of 0.9 was specified, and no orphans 


were generated in the overlap process for the baseline and 
refined-blade grids. The overlap process for the coarse grid 
generated a small number of orphans (on the order of a few 
dozen) during part of the rotor revolution. 

Computational Results 

Most results presented here are obtained using a standard 
set of input parameters for CFD modeling. This stan- 
dard or “best practices” set has emerged over the course of 
numerous rotorcraft simulations using the FUN3D solver. 
The standard computation utilizes a mixed-element mesh 
with prisms in the boundary layers near solid surfaces, and 
tetrahedra elsewhere. Inviscid fluxes are evaluated with a 
second-order (in space) implementation of Roe’s scheme 
(Ref. 18). The standard turbulence model is the Spalart- 
Allmaras turbulence model (Ref. 19). The turbulence model 
is applied to both the near-body and off-body regions. For 
the present results, no rotational corrections to the turbu- 
lence model are utilized. Time advancement is carried out 
with a modified second-order backward-difference dual- 
time scheme which results in reduced leading-order trunca- 
tion error compared to the classical second-order backward- 
difference scheme (Refs. 20, 21). The standard time step 
corresponds to a 1° change in blade azimuth per step. Time 
advancement is performed under the guidance of a ‘tempo- 
ral error controller’ (Refs. 11,20) in which dual-time subit- 
erations are carried out within each time step until either the 
L 2 norm of the nonlinear residual is reduced to a specified 
fraction of the estimated temporal error, or, a specified max- 
imum number of subiterations have been performed. For 
the results presented here, the error fraction ranges from 
0.005 to 0.01, and the maximum number of subiterations 
allowed ranges from 40 to 60. 

To assess the sensitivity of the computed results to sev- 
eral important CFD “parameters”, variations from the stan- 
dard set described above are examined. The effect of cell 
topology is considered by using both tetrahedral and mixed- 
element grids. The effect of mesh refinement is considered 
by performing isolated-rotor simulations using the coars- 



ened, baseline, and refined-blade grids described above, 
with resolutions ranging from 6.9 to 37.4 million nodes. 
The influence of the fuselage on the rotor airloads is ex- 
amined by replacing the empty-box background grid of the 
baseline isolated-rotor mesh with one containing the fuse- 
lage, leaving the rotor grids unchanged. Sensitivity to time 
step is examined using a step which corresponds to a 0.5° 
change in blade azimuth per step. Two additional turbu- 
lence models are considered: the Shear Stress Transport 
(SST) model of Menter (Ref. 22) and a hybrid model that 
combines the SST model near walls with a Large Eddy Sim- 
ulation (LES) model of the subgrid-scale turbulent kinetic 
energy (Ref. 23). 

The two flight conditions considered in this paper are de- 
noted as counters 8534 and 9020, and represent high-speed 
and high-thrust conditions, respectively. Table 3 summa- 
rizes the flight conditions and trim targets used for the sim- 
ulations. Within CAMRAD II, the thrust target is specified 
via the solidity-weighted thrust coefficient while the pitch 
and rolling moment targets are specified directly. The ro- 
tor is tilted relative to the fuselage by —3° (forward tilt, see 
Fig. 3). This tilt is accounted for in the angle of attack used 
in the simulation when running isolated rotor analyses. 

Figure 4 illustrates the radial stations at which flight- 
test data were collected on one of the rotor blades. Here, 
airloads are presented at a subset of these measurement 
stations, with comparison of flight data and computed re- 
sults for normal force, pitching moment, and chord force 
vs. blade azimuth at selected stations. The zero position for 
blade azimuth is defined as that for which the blade is nomi- 
nally aligned with the fuselage axis, pointing aft. Computa- 
tional results for forces and moments include both pressure 
and skin-friction contributions, whereas flight data only in- 
clude pressure contributions. The plots at each selected sta- 
tion are presented in a triptych format, which makes it easy 
to see the simultaneous effects of an ‘event’, such as dy- 
namic stall, on all three quantities. Sectional airloads at 
each radial station are presented in nondimensional form, 
with freestream speed of sound used as the reference veloc- 
ity - e.g. M 2 C n . Note that means are removed from the sec- 
tional airloads plots. It is known that there are discrepancies 
in some of the flight data due to faulty pressure taps. Plot- 
ting data with the means removed masks some of these dis- 
crepancies and facilitates assessment of the computational 
trends relative to flight. However, mean airloads (averaged 
around the rotor disk) are also compared as a function of 
radius in order to assess the agreement of the absolute lev- 
els. The sectional airloads plots are used to illustrate the 
sensitivity of the computed results to changes in CFD pa- 
rameters. To keep the number of figures to a manageable 
level, most plots are structured to show the sensitivity of 
the airloads (means removed) to the particular parameter 
under consideration at two of the nine radial stations, along 
with the sensitivity of the mean airloads as a function of 
radius. The two stations chosen are either the two showing 


the most sensitivity, or one showing the most sensitivity and 
one showing a typical level of sensitivity. 

Counter 8534: High-Speed Condition 

Flight counter 8534 has become a canonical case for 
CFD/CSD validation, so that condition is considered first. 
This flight condition represents the highest level-flight 
speed in the UH-60A Airloads Database, as well as the 
flight condition with the highest vibration levels. The 
salient features for this flight condition are given in Table 
3. 

To examine the effects of mesh size, computations on 
three different mixed-element meshes, containing between 
6.9 million to 37.4 million nodes, are shown in Figure 5 
for an isolated rotor. The top section of Figure 5 shows 
the variation of airloads with azimuth angle for station 
r/R = 0.965, while the mid section shows the variation at 
r/R = 0.990 and the bottom section illustrates the mean 
(azimuthally averaged) variation with radius. The level of 
differences observed at station r/R = 0.965 are typical of 
differences seen at other stations. By far the largest differ- 
ences are observed near the tip at r/R = 0.99. Mesh refine- 
ment has virtually no effect on the mean airloads. Similarly, 
time-step refinement, shown in Figure 6, shows no effect 
on the rotor airloads as the step is reduced from the stan- 
dard value of 1 .0° to 0.5°. In a recent paper by Sitaraman et 
al. (Ref. 24), using feature-based off-body adaption, a sim- 
ilar lack of sensitivity of rotor airloads to mesh refinement 
was observed for the high-speed (C8534) flight condition. 

It should be noted that trim targets for rotor thrust and 
pitch/roll moments used in the CFD/CSD coupling process 
are indirectly generated from flight data using estimates 
of the vehicle gross weight, tail plane loads, and fuselage 
loads. They are not obtained by integration of the measured 
sectional loads. Thus differences in the mean sectional air- 
loads and moments between computation and flight do not 
indicate that the specified trim targets are not met. 

Figure 7 shows the comparison of data and computed 
results for an isolated rotor obtained using both a fully- 
tetrahedral mesh and a mixed-element mesh. Both meshes 
have identical point distributions, approximately 17.6 mil- 
lion nodes. With means removed, the variation of the com- 
puted airloads with blade azimuthal angle generally agree 
well with the flight data. Other stations show similar lev- 
els of agreement between flight data and computation. Very 
little difference is observed in the computed results on the 
two meshes for this flight condition. The largest differences 
are seen near the rotor tip, where the mixed-element mesh 
result matches the data slightly better than the tetrahedral 
mesh. In fact at r/R = 0.990 the results on the 17.6 million 
node tetrahedral mesh are nearly identical to the results on 
the coarsest mixed-element mesh shown in Figure 5. For 
all other stations not shown, the differences between the 



Counter 

a (deg) 


M tip 

Ct / G 

My (ft-lb) 

M x (ft-lb) 

8534 

-4.31 

0.364 

0.642 

0.081 

-4169 

-6042 

9020 

2.49 

0.245 

0.669 

0.118 

-2176 

-616 


Table 3. Flight conditions and trim targets; sign conventions: a + nose up, My + nose up, Mx + right side down. 


tetrahedral mesh results and the mixed-element mesh re- 
sults are very small indeed, comparable to the differences 
observed at r/R = 0.865. Considering the mean airloads 
(averaged around the rotor disk) as a function of radius, the 
normal force is reasonably well predicted, the pitching mo- 
ment less so, and the chord force is poorly predicted past 
r/R = 0.40. The flight data for both pitching moment and 
chord force exhibit a significant amount of variation from 
station to station past r/R = 0.40. The computational re- 
sults exhibit much less station- to- station variation than the 
flight data, although the trends past r/R = 0.80 appear sim- 
ilar. The use of the mixed-element grid results in slightly 
better agreement in the mean chord force along the blade. 
The mean normal force and mean pitching moment are vir- 
tually indistinguishable between the two computations. 

The effect of the fuselage is shown in Figure 8 at the two 
innermost radial stations. Although not shown, the effect of 
the fuselage on the airloads for r/R > 0.5 is much less sig- 
nificant than at stations inboard. For the innermost stations, 
inclusion of the fuselage results in a general improvement 
in the correlation between computation and flight data over 
the forward part of the rotor disk, particularly for the chord 
force and pitching moment. No substantial change in the 
correlation with flight data over the aft portion of the rotor 
disk is obtained by including the fuselage; this area may 
be more strongly influenced by the hub (not included in the 
present computations). The mean airloads show virtually no 
sensitivity to the presence of the fuselage for the outboard 
sections of the blade, and only minimal differences for the 
inboard stations. 

Next, for an isolated rotor, using the baseline, mixed- 
element grid, variations of turbulence model (SA, SST and 
HRLES) are considered. Results are shown for stations 
r/R = 0.965 and r/R = 0.99 in Figure 9. Most surprisingly, 
none of the turbulence models, for this flight condition, pro- 
duced any noticeable change in the computed results. Al- 
though the mesh is fairly fine and relatively isotropic in the 
region of the rotor disk, it is an open question whether the 
grid is fine enough, and the time step small enough, to best 
exploit the the HRLES model. 

Table 4 summarizes the sensitivity of the blade control 
angles required for trim of the rotor to the various CFD pa- 
rameters discussed above. As expected from the lack of 
sensitivity of the airloads for this flight condition, the con- 
trol angles show very little change with parameter variation, 
a few tenths of a degree at most. In addition, Table 4 sum- 
marizes the sensitivity of the computed thrust and torque 


coefficients to CFD parameter variation. Taking the base- 
line, mixed-element mesh for the isolated rotor as the refer- 
ence result, over the range of parameters considered, thrust 
varies by no more than 0.3 percent, and torque varies by no 
more than 3.5 percent. 

Next the structural loads computed within the coupled 
CFD/CSD process are compared with those obtained from 
flight. The UH-60A database contains measured data for 
torsional moments, normal (flap) moments, and edgewise 
(chordwise) moments at several stations not coincident with 
the airloads stations. Figure 10 shows the comparison of 
computed and measured torsion moment at four radial sta- 
tions, while Figures 11 and 12 show normal and edge- 
wise bending moments, respectively. The results are shown 
for the baseline, mixed-element mesh for an isolated rotor, 
the corresponding tetrahedral mesh, and the mixed-element 
mesh with the fuselage, all with 1° time step and the SA 
turbulence model. Of the CFD parameter variations con- 
sidered, these three cases span the range of variation in the 
structural loads. As with the airloads, means have been re- 
moved. In Figure 10, the torsional moment at r/R = 0.0466 
has been obtained indirectly from the pitch-link load and the 
pitch-link load offset of 0.0224R; the same procedure was 
done for both flight and computation. In general, the 1/rev 
variation of torsional and normal-bending moments are rea- 
sonably well predicted by the computation. Higher harmon- 
ics are not well captured, especially for the flap-bending 
moment atr/R = 0.113. Overall the edgewise bending mo- 
ments are rather poorly predicted, even for the lowest har- 
monic. As suggested by the airloads predictions, the struc- 
tural loads show the most sensitivity to the inclusion of the 
fuselage, although the influence of the fuselage on the re- 
sults is small. 

Counter 9020: High-Thrust Condition 

Flight counter 9020 represents a high-thrust condition at a 
moderate advance ratio; the conditions for this counter are 
summarized in Table 3. The high-thrust level leads to dy- 
namic stall for several radial stations on the retreating side. 

Mesh refinement was also carried out for counter 9020. 
For this counter, unlike counter 8534, acceptable iterative 
convergence during each time step was typically not ob- 
tained using a step of 1° on the refined blade grid. For 
that mesh, with a step of 1°, convergence at nearly all time 
steps would “hang”, and fail to reach the specified tempo- 
ral error tolerance (0.01 times the residual norm); because 



Mesh 

Time Step 

Turb. Model 

0o 

01 c 

01s 

Cj 

Cq 

CO-mix-iso 

1.0° 

SA 

13.51 

2.47 

-8.99 

0.00688 

0.000703 

BL-mix-iso 

1.0° 

SA 

13.40 

2.42 

-8.94 

0.00688 

0.000682 

RB -mix-iso 

1.0° 

SA 

13.37 

2.42 

-8.95 

0.00688 

0.000675 

BL-tet-iso 

1.0° 

SA 

13.44 

2.49 

-8.97 

0.00688 

0.000706 

BL-mix-iso 

0.5° 

SA 

13.42 

2.42 

-8.97 

0.00688 

0.000683 

BL-mix-iso 

1.0° 

SST 

13.42 

2.43 

-8.95 

0.00688 

0.000681 

BL-mix-iso 

1.0° 

HRLES 

13.41 

2.43 

-8.96 

0.00688 

0.000685 

BL-mix-fuse 

1.0° 

SA 

13.51 

2.71 

-9.30 

0.00690 

0.000673 


Table 4. Sensitivity of control settings (measured at pitch hinge), thrust, and torque with respect to CFD parameter 
variation, Counter 8534. CO: 6.9 million nodes (coarse), BL: 17.6 million nodes (baseline), RB: 37.4 million nodes 
(refined blade). 


the convergence stalled, further subiterations did not lead 
to additional reduction in the residual. Reducing the time 
step to 0.5° rectified these convergence issues. As a result, 
mesh refinement was performed using a time step of 0.5°, 
to maintain an “apple-to-apples” comparison. 

Figure 13 shows the effect of mesh refinement at r/R = 
0.865 and r/R = 0.920 for an isolated rotor, along with the 
corresponding flight data. Two dynamic stall events are ob- 
served at these stations. Near T 1 = 270°, the lift drops pre- 
cipitously, accompanied by a large negative pitching mo- 
ment and reduction in chord force. After the first event, the 
lift builds again, until a second stall occurs near T 1 = 330°. 
In general, the phase of these stall events is in good agree- 
ment with the flight data, although the magnitudes (particu- 
larly pitching moment) differ from those observed in flight. 
Overall the effect of mesh refinement for counter 9020 is 
small, although some local effects are significant. Figure 
13 shows that while the first stall at r/R = 0.865 is un- 
changed when using the refined blade grid, the magnitude 
of the pitching moment drop is adversely impacted relative 
to the data for the second stall. At r/R = 0.920, the already 
delayed (relative to flight) second stall is delayed further, 
although the magnitude of the change in pitching moment 
when stall does occur is relatively unchanged. These two 
stations reflect the largest changes for this condition with 
mesh refinement; mean values are virtually unchanged. The 
results from Sitaraman et al. (Ref. 24) using feature-based 
off-body adaption also indicate a similar sensitivity of ro- 
tor airloads to mesh refinement for a related high-thrust 
(C9017) flight condition. 

Figure 14 shows that for the baseline mesh, temporal 
refinement from 1° per time step to 0.5° per time step pro- 
duced very little change in the computed results. The pri- 
mary effect of the smaller time step was to reduce some 
of the high-frequency oscillations observed after the dy- 
namic stall events in the computed results, particularly for 
the pitching moment. This is most evident at station r/R = 
0.920 in Figure 14. There is no significant effect on the 
mean airloads. 

Figure 15 shows computed results, for the isolated ro- 


tor, at stations r/R = 0.865, and r/R = 0.920 for the base- 
line mesh size using both mixed elements and pure tetra- 
hedra. In this high-thrust case, more differences appear 
between the fully-tetrahedral mesh and the mixed-element 
mesh than were observed in the high-speed case. Focusing 
on the pitching moment at r/R = 0.865, both grid topolo- 
gies capture the two stall events reasonably well. The tetra- 
hedral mesh leads to the first stall event occurring slightly 
early, while with the mixed-element mesh the onset of stall 
matches the flight data. Both meshes show a larger drop 
in pitching moment than the data, with the mixed-element 
mesh having a slightly larger reduction in pitching moment 
during the stall. For the second dynamic stall, the trend with 
mesh topology is reversed - the tetrahedral mesh exhibits a 
larger reduction in pitching moment. Considering station 
r/R = 0.920, the mixed-element mesh leads to a first stall 
that is in very good agreement with the data; the tetrahe- 
dral mesh again results in an early stall, with a significantly 
larger drop in pitching moment. As for station r/R = 0.865, 
the tetrahedral mesh results in the second stall occurring vir- 
tually at the same azimuthal location as flight, with reason- 
ably good agreement with flight for the magnitude of the 
drop. On the other hand, the mixed-element mesh results in 
the second stall occurring approximately 12° too late, with 
a noticeably smaller reduction in moment. 

Turning to the variation of the mean airloads with radius, 
the two meshes lead to nearly identical results for normal 
force, with slightly larger differences in pitching moment 
and chord force. The peak normal force, which occurs near 
r/R = 0.865 in the data, is underpredicted in the compu- 
tations. It would appear that there are some anomalies in 
the flight data near r/R = 0.400, since the normal force at 
that station is significantly different than neighboring radial 
locations, and the chord force there is off the scale of the 
plot. Interestingly, the flight data for pitching moment at 
r/R = 0.400 does not appear to be out of line with nearby 
radial locations. 

Next, the effects of including the fuselage are consid- 
ered. Compared with the high-thrust case, including the 
fuselage has a generally smaller effect on the results, but 



including the fuselage does not consistently improve the 
comparison with data. For example, Figure 16 shows the 
effect of the fuselage on the airloads at r/R = 0.400 and 
r/R = 0.675. At r/R = 0.400, inclusion of the fuselage 
improves the correlation with the flight data for the nor- 
mal force and chord force. While inclusion of the fuse- 
lage changes the pitching moment, the addition cannot be 
said to affect the overall level of agreement with flight. At 
r/R = 0.675, including the fuselage results in a markedly 
reduced drop in pitching moment near i jf = 240°, while the 
overall comparison with normal force and chord is neither 
better nor worse. The computed mean airloads for this flight 
condition are virtually unchanged when the fuselage is in- 
cluded. 

Because mesh refinement for flight counter 9020 showed 
some sensitivity to mesh size, turbulence model variations 
were conducted on the refined blade grid. Furthermore, 
because the refined blade grid exhibited poor convergence 
with the larger time step, a time step corresponding to 
0.5° azimuth change per step was used for the turbulence 
model variations. Figure 17 shows the comparison of the 
SA, SST, and HRLES models at stations r/R = 0.865 and 
r/R = 0.920. These radial stations show the most varia- 
tion between models over all the nine flight-data stations. 
At r/R = 0.865, all three models pick up the proper phase 
for the first stall. The SST model does a very good job at 
predicting the magnitude of the pitching moment through 
the first stall, while both SA and HRLES predict a larger 
negative pitching moment than indicated by the flight data. 
At r/R = 0.865 all three models predict the second stall to 
some degree, although once again it is the SST model that 
does the best job with the level of pitching moment change. 
At r/R = 0.920, only the SST model shows the slightest 
hint of the second stall event. The mean airloads are virtu- 
ally unaffected by the various turbulence models. 

The sensitivity to CFD parameters of the trim control 
angles and the rotor thrust and torque are listed in Table 
5. Compared to the high-speed case, more variation in the 
results for the torque coefficient are observed when the pa- 
rameters are varied. Taking the baseline, mixed-element 
mesh for the isolated rotor as the reference result, over the 
range of parameters considered, thrust varies by no more 
than 0.3 percent, whereas torque varies by as much as 13.5 
percent. 

Finally, Figure 18 shows the comparison of computed 
and measured torsion moment at four radial stations, while 
Figures 19 and 20 show the normal and edgewise bending 
moments, respectively. The results are shown for the base- 
line, mixed-element mesh for the isolated rotor, the corre- 
sponding tetrahedral mesh, and the tetrahedral mesh with 
the fuselage, all with 1° time step and the SA turbulence 
model. As for the high-speed case, overall variation is rea- 
sonably well predicted for the torsion and normal bending 
moments, with higher-harmonics being less well predicted, 
especially for normal bending at r/R = 0.113. Edgewise 


bending shows little agreement with data, although for this 
flight condition, clipping is evident in the flight data. As 
expected from the airloads results, little sensitivity to CFD 
parameter variation is evident in the structural loads. 

Summary 

A fully unstructured RANS solver has been used to com- 
pute rotor airloads for the UH-60A helicopter, for the pur- 
pose of validating the solver against the flight database. 
Two flight conditions, corresponding to high-speed and 
high-thrust conditions were examined. Parametric varia- 
tions of the computed results to mesh size, time step, and 
turbulence model were considered in the CFD model; the 
structural model was held fixed in all cases. Remarkably, 
the high-speed case exhibited virtually no sensitivity to any 
of the variations in computational inputs, except for the in- 
clusion of the fuselage. Although the unsteady (i.e. mean 
removed) variation of the computational results showed rea- 
sonably good agreement with flight data, there were dif- 
ferences between computation and flight that were not im- 
pacted by variation of the computational input. Further- 
more, the agreement between computation and flight for 
the mean loads was not nearly as good as for the unsteady 
loads. Likewise, the agreement between computation and 
flight for the mean loads was not significantly influenced by 
computational variation. There are a number of results from 
CFD/CSD computations for Counter 8534 in the literature, 
for example fully structured-grid results of Reference 1 and 
the hybrid unstructured/Cartesian results of Reference 24. 
The fully unstructured results presented in this paper are 
virtually identical to results in the cited references. While 
the high-speed case showed no sensitivity to CFD-model 
variations, the high-thrust case did show some sensitivity, 
the largest being with respect to the second of two dynamic 
stall events on the retreating side. 

The relatively small changes in computational results 
with changes in important CFD modeling parameters (mesh 
size, turbulence model, and time step) may suggest that the 
observed differences with flight data are not wholly due 
to errors in CFD modeling, and that perhaps the structural 
modeling needs to be examined. Of course there are still ar- 
eas lacking in CFD modeling. While a reasonably detailed 
model of the fuselage was included in the computations, 
the rotor hub and associated linkages, were not modeled, 
nor was the tail rotor. These missing components may well 
influence the flow over the rear of the rotor disk. However, 
they are unlikely to account for discrepancies between com- 
putation and flight elsewhere on the rotor disk. 
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Mesh 

Time Step 

Turb. Model 

0o 

Ole 

01* 

Cj 

Cq 

CO-mix-iso 

0.5° 

SA 
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3.01 
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0.000540 

BL-mix-iso 

0.5° 

SA 
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BL-tet-iso 

1.0° 

SA 
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RB -mix-iso 
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SA 
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RB -mix-iso 
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SST 
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RB -mix-iso 

0.5° 

HRLES 
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1.0° 

SA 
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3.07 
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Table 5. Sensitivity of control settings (measured at pitch hinge), thrust, and torque with respect to CFD parameter 
variation, Counter 9020. CO: 6.9 million nodes, BL: 17.6 million nodes, RB: 37.4 million nodes 
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(a) Coarse surface grid 


(b) Baseline surface grid 


(c) Refined-blade surface grid 





(d) Coarse volume slice 


(e) Baseline volume slice 


(f) Refined-blade volume slice 


Fig. 1. Comparison of UH-60 blade surface grids and volume slices at r/R = 0 . 865 . 



Fig. 2. Surface grids corresponding to the 19 million 
node baseline volume mesh. 


Fig. 3. Slice through the baseline volume grid systems, 
with and without fuselage (19 million and 18 million 
grid points, respectively). 



r/R=0.965 


Fig. 4. Blade radial stations for which flight airloads are available. 
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Fig. 5. Effect of mesh refinement on computed airloads (isolated rotor) at r/R = 0.965, and r/R = 0.990, together 
with azimuthal averages, Counter 8534. CO: 6.9 million nodes, BL: 17.6 million nodes, RB: 37.4 million nodes; all 
mixed-element meshes. 
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Fig. 6. Effect of time step refinement on computed airloads (isolated rotor) at r/R = 0.965, and r/R = 0.990, together 
with azimuthal averages, Counter 8534. 
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Fig. 7. Comparison of computed (isolated rotor) and measured airloads at r/R = 0.865, r/R=0.990, together with 
azimuthal averages, Counter 8534. Computations performed on mixed-element and tetrahedral meshes. 
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Fig. 8. Effect of fuselage on computed airloads at r/R = 0.225, and r/R = 0.400, together with azimuthal averages, 
Counter 8534. Computations performed on mixed-element meshes. 
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Fig. 9. Effect of turbulence model variation on computed airloads (isolated rotor) at r/R = 0.965, and r/R = 0.990, 
together with azimuthal averages, Counter 8534. Computations performed on the 17.6 million node mixed-element 
mesh, with a time step corresponding to 1.0°. 
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Fig. 10. Effect of selected CFD mesh variations on the computed torsional moment (means removed), Counter 8534. 




Fig. 11. Effect of selected CFD mesh variations on the computed normal bending moment (means removed), 
Counter 8534. 
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Fig. 12. Effect of selected CFD mesh variations on the computed edgewise bending moment (means removed), 
Counter 8534. 
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Fig. 13. Effect of mesh refinement on computed airloads (isolated rotor) at r/R = 0.865, and r/R = 0.920, together 
with azimuthal averages, Counter 9020. CO: 6.9 million nodes, BL: 17.6 million nodes, RB: 37.4 million nodes). 
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Fig. 14. Effect of time step refinement on computed airloads (isolated rotor) at r/R = 0.865, and r/R = 0.920, together 
with azimuthal averages, Counter 9020. Computations performed on the 37.6 million node mesh 
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Fig. 15. Comparison of computed (isolated rotor) and measured airloads at r/R = 0.865, and r/R = 0.920, together 
with azimuthal averages, Counter 9020. Computations performed on mixed-element and tetrahedral meshes. 
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Fig. 16. Effect of fuselage on computed airloads at r/R = 0.400, and r/R = 0.675, together with azimuthal averages, 
Counter 9020. 
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Fig. 17. Effect of turbulence model on computed airloads (isolated rotor) at r/R = 0.865, and r/R = 0.920, together 
with azimuthal averages, Counter 9020. Computations performed on the 37.6 million node mixed-element mesh, 
with a time step corresponding to 0.5°. 




Fig. 18. Effect of selected CFD mesh variations on the computed torsional moment (means removed), Counter 9020. 




Fig. 19. Effect of selected CFD mesh variations on the computed normal bending moment (means removed), 
Counter 9020. 
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Fig. 20. Effect of selected CFD mesh variations on the computed edgewise bending moment (means removed), 
Counter 9020. 


